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PG 1018-047: the longest period subdwarf B binary 
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ABSTRACT 

About 50% of all known hot subdwarf B stars (sdBs) reside in close (short period) binar- 
ies, for which co mmon envelope ejection is the most likely formation mechanism. However, 
iHan et all (|2()03|) predict that the majority of sdBs should form through stable mass transfer 
leading to long period binaries. Determining orbital periods for these systems is challenging 
and while the orbital periods of 100 short period systems have been measured, there are no 
periods measured above 30 days. As part of a large program to characterise the orbital periods 
of subdwarf B binaries and their formation history, we have found that PG 1018-047 has an 
orbital period of 760 ± 6 days, easily making it the longest period ever detected for a subd- 
warf B binary. Exploiting the Balmer lines of the subdwarf primary and the narrow absorption 
lines of the companion present in the spectra, we derive the radial velocity amplitudes of both 
stars, and estimate the mass ratio Mms / M^dB = 1-6 ± 0.2. From the combination of visual 
and infrared photometry, the spectral type of the companion star is determined to be mid K. 

Key words: binaries: close - binaries: spectroscopic - subdwarfs - stars: evolution - stars: 
individual: PG 1018-047 



1 INTRODUCTION 

Subdwarf B stars (sdBs) are core-helium burning stars with thin 
hydrogen envelopes. They are situated between the main sequence 
and the white dwarf cooling track at the blueward extension of the 
horizontal branch, the so-called Extreme or Extended Horizontal 
Branch (Heberetal. 1984; Heber 1986; Safferetal. 1994). Subd- 
warf B stars have colours and spectral characteristics correspond- 
ing to those of a B star, but the Balmer lines are abnormally broad 
for the colour compared with population I main-sequence B stars 
due to their high surface gravities (logg ~ 5.0 — 6.0). Subd- 
warf B stars have a typical mass of 0.5 Mq i lHeberetalJ[T984h 
and can be found in all Galactic populations. They are thought 
to b e the dominant source for the UV-uptum in early type galax- 
ies jperguson et al. ^991; Bro wn et alj 2000l). A fra ction of sdBs 
pulsate jCharpinet et al.lil996i ; 'Kilkenny et al1ll997l) . giving great 
opportunities to derive fundamental paramete rs (e.g. the stellar 
mass) and study their internal structure in d etail dOreen et al .120031 : 
iFontaine et ai](2008l : l0stensenll2OO^ I2OI0I) . Subdwarf B stars are 
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also suggested to be ver y useful as age indi cators using evolution- 
ary population synthesis (iBrown et al.l 19971) . or as distance indicat- 
ors jKilkennv etd]| 19991) 7 Since a large fraction of sdBs are mem- 



bers of binary systems (^ Maxted et al. 1200 ll) and because they are 
intrinsically bright and ubiquitous, they are therefore an ideal pop- 
ulation in which to study binary star evolution. For a comprehensi ve 
review on hot subdwarf stars we refer the reader to lHebed l l2009h . 

iHan et all |2003) describe in detail the formation and evolu- 
tion of sdBs by using binary population synthesis models. They 
find that sdBs form via five main evolutionary channels: the first 
and second common envelope channels, the first and second stable 
Roche lobe overflow channels and the helium white dwarf merger 
channel. This last channel is the only one that results in the form- 
ation of single sdBs. They find that the contribution of the second 
Roche lobe overflow channel is not significant, leaving only three 
channels to form sdB binaries. Each of these three binary forma- 
tion channels predicts a different orbital period distribution for the 
population of sdBs. The binaries formed through the first common 
envelope channel should display orbital periods between 0.5 and 
^40 days and the companions to the sdBs will be main sequence 
stars. Binaries formed via the second common envelope channel 
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are expected to have white dwarf companions and their range of or- 
bital periods will be wider, extending further into the short periods 
but not to long periods. Note as well that these common env elope 
phases are not very well understood. iNelemans et"^ feOOOh have 
concluded from the observed double white dwarf population that 
its outcome may not always be a strong reduction of the orbital sep- 
aration. Finally, sdB binaries formed through the first stable Roche 
lobe overflow channel will have main sequence companions and 
will display orbital pe riods between 0.5 and ~2000 days. 

iHan et al ] l l2003h conclude that their set 2 of simulations is 
the model th at best describ es the observed sample of short period 
sdB binaries ( Moral es-Rueda et al.il2003l) . In this particular model 
(and also in 9 out of the 12 models they describe) the majority of 
sdB binaries, between 60 and 70 percent of the total, are formed 
via the first stable Roche lobe overflow channel. At the same time, 
this is the channel most affected by observational selection effects 
decreasing the number of observable sdB binaries formed through 
this channel. These observational effects are primarily that sdBs 
with companions that are brighter than the sdB itself will not be 
identified as sdBs at all. The second effect has to do with obser- 
vational limitations: it is easier to detect radial velocity variations 
from a short orbital period system than from a long one, as these 
are smaller and take longer to determine in the second case. 

Despite extensive observational work, not a single system has 
been found in this long-period regime (first stable RLOF channel), 
wherea s at present ~100 sd B binaries with short peri ods are con- 
firmed jGeier et al.l201ll and lCopperwheat et alj201]]) . The orbital 
periods are mostly below 1 day, with a median period of 0.61 days. 
However, in this work we will report on PG 1018-047, the first truly 
long-period subdwarf B binary. Whi le this may be the pr oduct of 
binary evolution as sugg ested by the lHan et all ( l2002l.[20b3) mod- 
els, we will find that it could also be the remnant of a hierarchical 
triple, as recently outlined bv'Claus en & Wad3 ( 1201 ih . 

Our target, PG 1018-047 has an apparent visual 
(Stromgren) magnitude niy — 13.32 and was discovered as 
an ultravi olet-excess stella r obje ct in the Palomar-Green Survey 
dGreen, Sc hmidt & Lieberj Il986l) . It was subsequently observed 
by Maxt ed et al.l ( 120011) to check for radial velocity variations. 
Although weak spectral f eatures from a late-t ype companion are 
visible in the spectrum, iMaxted et alj ( 1200 ih did not find any 
significant radial velocity shifts using their variability criteria. This 
led to the conclusion that PG 1018-047 is probably not a binary 
with a short orbital period. The presence of the companion in the 
spectrum prompted continued follow-up of the system in order to 
determine how such a binary could have formed in the first place. 
We present the results after more than a decade of monitoring. 



2 OBSERVATIONS AND REDUCTION 

We have observed PG 1018-047 spectroscopically with several dif- 
ferent instrument setups over a period of ten years. In Table [T] we 
summarize the observing dates, the setup used in each case, the 
wavelength range covered and the number of spectra obtained dur- 
ing each epoch. The data were obtained using the Isaac Newton 
Telescope (INT), William Herschel Telescope (whtQ and Nordic 
Optical Telescope (NOT) on the Island of La Palma, the Radcliffe 

^ Both the INT and WHT belong to the Isaac Newton Group of Telescopes 
(ESfG). 



telescope at the South African Astronomical Observatory (SAAO) 
and the Hobby-Eberly Telescope (HET) located at the McDonald 
Observatory in Texas. The different instrimient setups were as fol- 
lows: 

For the INT-Red spectra the intermediate dispersion spectro- 
graph (IDS) was used. It is a long-slit spectrograph mounted on 
the Cassegrain Focal Station of the INT. The 500 mm camera to- 
gether with the high resolution R1200R grating and a windowed 
Tek5 CCD centered in A = 6560 A covered the Ha region. A 
I arcsec slit was used. 

INT-Blue: The INT with the IDS, equipped with the 235 cam- 
era, the R1200B grating and a windowed EEVIO CCD was used to 
obtain these blue spectra. The 2002-2007, 2008-2009 spectra were 
centered on respectively A = 4348 A and A — 4505 A, cover- 
ing as many B aimer lines to the blue as possible, including Hb at 
4861.327 A. For all exposures a 1 arcsec slit was used. 

The SAAO data were obtained using the Radcliffe 1.9 m tele- 
scope together with the grating spectrograph plus the SITe back- 
illuminated CCD. Grating 4, with 1200 grooves per millimeter was 
used to obtain spectra covering and Hp with a dispersion of 
0.5 A/pix and a resolution of 1 A at 4600 A. The slit width varied 
from 1.2 to 1.5 arcsec depending on the seeing. 

WHT-Red: The WHT was equipped with the double arm In- 
termediate dispersion Spectrograph and Imaging System (ISIS). 
The R1200R grating and the Red+ CCD were used to obtain the red 
spectra centered on A = 6560 A (2007 data) and on A = 6521 A 
(2009 data). A slit width of 1.2 arcsec was used for the 2007 obser- 
vations and a 1 arcsec slit for the 2009 ones. 

The setup WHT-Blue denotes WHT data obtained using the 
ISIS spectrograph with the R600B grating and the blue EEVIO. 
The grating was centered on A = 4388 A with a 1 arcsec slit (2006 
observations), on A = 4500 A with a 1.5 arcsec slit (2007 data), 
A = 4339 A with a slit width of 0.62 arcsec (2008) and on A = 
4349 A with a 1.04 arcsec slit during the 2009 observations. 

The NOT data were taken using the 2.56 m Nordic Optical 
Telescope. The FIES (Fibre-fed Echelle Spectrograph) highest res- 
olution fiber (the 1.3 arcsec fibre offering a spectral resolution of R 
= 67000) covered the entire spectral range 3700 - 7300 A without 
gaps in a single fixed setting. 

The final setup, HET, refers to data from the bench-mounted 
echelle fibre-fed High Resolution Spectrograph, mounted on the 
9.2m Hobby-Eberly Telescope operated in its R = 15000 resolution 
mode. In the "2x3" on-chip binning mode that was used, the dis- 
persion was about 6 km/s per binned output pixel. A 2-arcsecond 
optical fibre was used for the stellar target, and two additional fibers 
were used to record the sky spectrum. A cross-dispersing grating 
with 600 grooves per millimeter was used, centering A ~ 5822 A 
at the boundary between the "blue" and "red" CCDs. The useful 
wavelength coverage extended from 4810 A to 6760 A. 

To reduce the spectra from the INT, WHT and SAAO, stand- 
ard Starlink routines were used. Flatfields were taken to correct for 
the pixel to pixel variations in the CCD and the bias correction 
was carried out by using the overscan region in each CCD frame. 
The ob j ects w ere extracted with the optimal extraction algorithm of 
lMarsri ( ll989l) . CuAr-l-CuNe arc spectra were taken before and after 
each target spectrum or after each set of two spectra at the target's 
position to calibrate these in wavelength. Fourth order polynomials 
were computed to fit the lines in the arcs and the solutions were 
used for the calibration of the corresponding spectra. 
The NOT data were reduced with the automatic data reduction soft- 
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Table 1. Journal of observations. Observers: P. F. L. Maxted (P. M.), T. Augustijn (T. A.), T. R. Marsh (T. M.), Luisa Morales-Rueda (L. M.), 
G. Nelemans (G. N.), C. Copperwheat (C. C), R.A. Wade (R. W.) and M. A. Stark (M. S.). 



Date 


Setu 
eup 


A region 


Spcctfti 


A4ean dispersion 
(A/pixcl) 


Ot)scrvGr(s) 


11 - 19/04/00 


INT-R 


Ha 


4 


0.39 


PM. 


08 - 13/03/01 


INT-R 


Ha 


11 


0.39 


PM. 


01 - 07/05/01 


INT-R 


Ha 


9 


0.39 


PM. 


26 - 30/03/02 


SAAO 


Blue 


6 


0.49 


T.M. 


25 - 27/04/02 


INT-B 


Blue 


6 


0.48 


T. A. & T. M. 


09 - 16/04/03 


INT-B 


Blue 


27 


0.48 


T.M. 


30/03 - 05/04/04 


SAAO 


Blue 


5 


0.49 


L. M. 


23 - 24/06/05 


SAAO 


Blue 


2 


0.50 


L. M. 


06/02/06 


WHT-R 


Ha 


2 


0.22 


G. N. 




WHT-B 


Blue 


2 


0.44 


G. N. 


09/03/07 


WHT-B 


Blue 


4 


0.44 


service 


27/03 - 07/04/07 


INT-B 


Blue 


1 1 


0.48 


T. M. 


29- 31/03/07 


WHT-R 


Ha 


5 


0.25 


G. N. 




WHT-B 


Blue 


5 


0.44 


G. N. 


21 - 22/03/08 


INT-B 


Blue 


4 


0.48 


C. C. 


01/05/08 


WHT-R 


Ha 


2 


0.49 


PM. 




WHT-B 


Blue 


2 


0.44 


PM. 


11/03/09 


INT-B 


Blue 


2 


0.48 


C. C. 


30/04/09 


WHT-R 


Ha 


4 


0.25 


T.M. 




WHT-B 


Blue 


4 


0.44 


T.M. 


03/04/10 


NOT 


Ha + Blue 


1 


0.03 


service 


06/12/07 - 23/03/10 


HET 


Ha + Blue 


7 


0.12 


R. W&M. S. 



be identified with Oil lines and possibly also Nil. We made a fit 
to the mean spectrum of PG 1018-047 using the LTE model grids 
of Heber et al. (2000), with explicit metals of solar composition 
and abundances depleted by 0.0, 0.5, 1.0, 1.5 and 2.0 dex relat- 
ive to solar. A reasonable fit is achieved with T^a = 30500 ± 200 
K, log g = 5.50 ±0.02, and with the N and O abundance 1/10 of 
the solar value. In Figure 1 , a synthetic spectrum with these para- 
meters broadened to match the resolution of the observed spec- 
trum is shown. The model spectrum contains helium at a fraction 
log N(He)/N(H) = -3.0, but even this is clearly too much. In order 
to make helium fit with the observed spectrum, the model must be 
depleted to log N(He)/N(H) < - 4. It is also clear that the abund- 
ances of the various elements are quite far from solar composition 
relative to each other, which is not unusual for the sdBs (Heber et 
al. 2000). Note that the K-star contributes some light also in the 
blue part of the spectrum, but insufficient to make any lines clearly 
visible in the spectrum. However, the contribution to the continuum 
might still be sufficient to affect the fitting procedure. In the H^ re- 
gion, however, metal lines from the K-star are clearly seen. Our 
high S/N mean spectrum has too low resolution to reliably infer the 
abundances of the individual components, and our high resolution 
spectra have insufficient S/N. 

3.2 The orbit of PG 1018-047 

3.2. 1 Radial velocity measurements 

The radial velocities (RVs) of the INT, WHT, SAAO and NOT 
ever, tne Jrlel line at 44 /z A, wnicn is typical ol sats stars is ao- spectra were determ i ned f ollowing the procedure described by 
sent. Instead the Si III triplet (4553, 4568, 4575 A) is the strongest iMorales-Rueda et"al] (2003*), i.e. least squares fitting of a line pro- 
feature. Several lines in the region between 4600 and 4700 A can file model. This line profile model was built up from three Gaus- 

sians per B aimer line with different widths and depths, but with 
a common central wavelength position which varies between the 
2 Developed by Eric Stempels (http://www.not.iac.es/iiis tru-rnents/fies/fiestool/FIE^eotianT tie parameters of the Gaussians were optimized by com- 



ware package FIEStooJf], which makes use of the IRAF and Nu- 
mArray packages via a Python interface. After preprocessing, the 
raw frames were debiased and subsequently divided by a 2D nor- 
malized Hatfield, correcting for the shape of each spectral order. 
Next the science spectra were extracted using the optimal extrac- 
tion algorithm from Home ( 1986) and corrected for the blaze shape. 
The wavelength calibration was done from a ThAr lamp spectrum 
taken right before the science data. For the HET spectra we used 
standard IRAF tasks, organized using"pipeline" scripts, to process 
the images and extract the spectra. Observations of PG 10 18-047 
were taken as pairs of 750-second exposures. Each pair was com- 
bined within IRAF, using the "crreject" option in the task "imcom- 
bine" to reduce cosmic-ray contamination; additional rejection of 
cosmic-ray artifacts was done later by hand. As a result of this ob- 
servational effort, we have a total of 125 spectra of PG 1018-047. 

3 RESULTS 

3. 1 The optical spectrum of PG 1018-047 

The normalised optical spectrum of PG 1018-047 is shown in Fig- 
ure[T] It is obtained by averaging over a number of the ING/SAAO 
spectra taken at the same orbital phase, calculated using our best 
orbital solution from section [T2l In order to normalise the spectra, 
we fitted third order polynomials to the regions free from absorp- 
tion lines and divided the spectra by these fits. The mean spectrum 
shows several metal lines in the region between H^ and H-, . How- 
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Figure 1. The mean spectram of PG 1018-047 for the region covering to Hi (left) and the region around Hq (light). Plotted below the 
mean spectrum (shifted down for clarity) is a model spectrum for an sdB star (see text for details). Lines from the K-dwarf are clearly seen in 
the red part, but the lines in the blue are from the sdB. 



Table 2. Narrow absorption lines pre sent in th e optical spectrum (Fig ure [T). The identification was done using the mid- and high-resolution 
spectral library of lMontes et al.l rl997l) and also lRalchenko et aDi201(ll . The lines with a * were included in the template to determine the radial 
velocities of the secondary (see section [3.2.3t . For reference, we have added the wavelengths for the Balmer lines we adopted in our analysis 
as well. 



A (A) Element A (A) Element A (A) Element A (A) Element 



4153.30 


On 


4596.17 


On 


6337.28* 


6609.05* 




4164.79 


Fe ni, blend 


4630.54 


Nil 


6359.45* 


6614.42* 




4189.80 


On 


4639.70 


On, blend 


6363.71* 


6637.09* 




4253.59 


Siii-hOii 


4642.26 


On, Nn 


6394.45* 


6644.39* 




4276.74 


On, blend 


4649.14 


On 


6400.43* 


6664.37* 




4414.91 


On 


4661.04 


On 


6408.96* 


6678.97* 




4416.58 


On 


4676.23 


On 


6412.37* 






4442.49 


On 


4700.31 


On, blend 


6422.05* 


4101.735 


Hs 


4552.65 


Sim 


4705.44 


On 


6431.44* 


4340.465 




4567.87 


Sim 


4710.04 


On 


6439.70* 


4861.327 




4574.78 


Si ni 






6450.51* 


6562.800 


Ha 


4590.57 


On 


6335.83* 




6456.49* 







paring the mo del to the normalized average spe ctrum over all ob- 
servations; see lMaxted. Marsh & Moraril ( l2000cl) for further details 
of this procedure. For the blue spectra we fit simultaneously for 
Hp, Hry and Hs, whereas for the red spectra only the Ha line can 
be fitted. The RV of the NOT spectrum was determined using a 
model containing Ha, Hp and H-y (see also Table |2j. The radial 
velocities from the 7 HET spectra are obtained from the Hp ab- 
sorption line using simple Gaussian fitting to the core of the line 
within the IRAF "splot" task. A list of the radial velocitities and 
the uncertainties measured is given in Table [3] 



3.2.2 Orbital parameters 

Once the radial velocities for all spectra were known, we use d 
the floating mean periodogram jCumming. Marcv & Butle r 1999), 
a generalizat i on of the we ll-known Lomb-Scargle periodogram 
( lLomb|[T976l : IScargldri982l) . to determine the most probable fre- 
quencies (periods) present in the data. The method consists in fit- 
ting the radial velocity data with a model composed of a sinusiod 



plus a constant of the form: 

v = 'y + Ksm[2'Kf{t-to)\, (1) 

with / the frequency and t the time of observation. This means 
we assume the binary system to have a circular orbit with semi- 
amplitude K and a systemic velocity 7. For each frequency 
(f=\IP) we perform least squares fitting of the data, solving 
for 7 and K sim ultaneously using singular value decomposition 
jPress et al.ll2002h . In this way we can obtain the statistic of the 
model as a function of frequency, or in other words the periodo- 
gram. 

In our initial period determination, the values turned out 
larger than expected given the number of data points, indicating 
that there must have been an extra unaccounted source of uncer- 
tainty, most likely due to systemic e ffects, intrinsic variability o f 
the star or slit-filling errors (see also iMorales-Rueda et ai] |2003h . 
Such errors are unlikely to be correlated with either the orbit or the 
statistical errors we have estimated. To allow for this we compute 
the level of systematic uncertainty (a) per telescope that when ad- 
ded in quadrature to our error estimates gives a reduced x^ ~ 1 
for each of the INT, WHT subsets, telescope by telescope, relat- 
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Table 3. The 125 radial velocity measurements for PG 1018-047 with their formal eiTors from the least squares fitting routine. 



HJD 


RV 


HJD 


RV 


HJD 


RV 


-2450000 


(km/s) 


-2450000 


(km/s) 


-2450000 


(km/s) 


1646.47135 


31.3 ± 4.4 


2741.52795 


49.8 ± 3.0 


4189.54710 


52.1 ± 1.4 


1646.47650 


20.5 ± 4.3 


2741.53909 


53.3 ± 3.1 


4189.54712 


42.9 ± 1.0 


1654.46074 


24.7 ± 3.8 


2741.54857 


48.8 ± 3.4 


4190.40980 


51.9 ± 1.0 


1654.46780 


32.1 ± 3.7 


2741.55807 


48.2 ± 4.3 


4190.40982 


46.0 ± 0.7 


1977.48220 


49.6 ± 5.3 


2743.42040 


50.4 ± 4.5 


4190.42420 


53.5 ± 1.0 


1978.62533 


55.2 ± 3.8 


2743.43451 


60.3 ± 4.5 


4190.42422 


50.7 ± 0.7 


1979.54044 


45.0 ± 3.2 


2743.52128 


52.2 ± 6.3 


4191.59677 


52.1 ± 1.4 


1979.59958 


65.5 ± 3.3 


2744.35895 


54.3 ± 3.5 


4191.59677 


44.1 ± 1.0 


1979.60655 


54.0 ± 3.3 


2744.37308 


49.7 ± 2.8 


4193.37192 


51.4 ±7.8 


1982.51489 


52.2 ± 3.9 


2744.39186 


49.4 ± 2.8 


4195.35718 


41.6 ± 1.9 


1982.52185 


54.6 ± 4.1 


2744.40598 


51.4 ± 2.7 


4195.37121 


36.4 ± 1.7 


1982.55421 


57.8 ±4.0 


2744.42529 


53.9 ±2.6 


4195.61447* 


35.5 ± 2.5 


1982.56118 


52.6 ± 3.9 


2744.43941 


52.2 ± 2.6 


4195.62457* 


38.1 ±4.5 


1982.60059 


57.5 ± 3.9 


2745.37982 


48.8 ± 3.5 


4197.51456 


43.8 ± 1.8 


1982.60755 


53.2 ± 4.0 


2745.39394 


45.3 ± 3.1 


4197.52858 


43.6 ± 1.9 


2031.36769 


62.1 ± 8.3 


2745.40987 


51.4 ± 2.8 


4441.00811 


39.7 ± 2.0 


2031.37466 


41.4 ± 8.4 


2745.42399 


51.9 ± 2.7 


4469.93158 


36.0 ± 2.0 


2031.38641 


49.5 ± 4.4 


2746.44131 


54.6 ± 2.5 


4502.85094 


33.5 ± 2.0 


2032.48243 


54.9 ± 3.0 


2746.45544 


51.0 ±2.5 


4547.37872 


33.0 ± 1.7 


2032.49633 


53.6 ± 3.5 


2746.47412 


55.4 ± 2.2 


4547.40665 


31.1 ± 1.6 


2033.39672 


53.8 ± 4.1 


2746.48824 


44.8 ± 2.3 


4547.62144 


29.6 ± 2.0 


2033.40369 


53.8 ± 4.3 


3095.38227 


26.8 ± 6.7 


4547.64844 


32.4 ± 2.2 


2037.45584 


54.1 ± 3.8 


3095.40336 


19.0 ± 6.6 


4550.71405 


26.7 ± 2.0 


2037.46282 


50.9 ± 3.5 


3098.33704 


12.9 ±9.8 


4562.68846 


29.6 ± 2.0 


2360.36683* 


-33.8 ± 25.0 


3101.31890 


18.8 ± 9.2 


4588.35123 


29.7 ± 2.5 


2360.37745* 


-4.6 ± 22.9 


3101.32956 


19.7 ± 8.7 


4588.35125* 


34.2 ± 2.6 


2362.36641* 


17.5 ± 19.3 


3545.21953 


51.4 ± 7.2 


4588.35646* 


22.8 ± 2.6 


2362.38049* 


5.4 ± 20.3 


3546.20756 


52.6 ± 5.8 


4588.35648 


29.5 ± 2.4 


2364.39275 


23.1 ± 6.0 


3772.55466 


35.6 ± 0.8 


4902.42320 


34.8 ± 4.5 


2364.40338 


17.9 ± 6.7 


3772.55467 


26.7 ± 1.2 


4902.44418 


36.6 ± 4.9 


2390.51113 


34.6 ± 17.1 


3772.68564* 


15.3 ±0.9 


4952.38360 


55.2 ±2.8 


2390.52179 


32.8 ± 9.0 


3772.68565 


34.6 ± 1.4 


4952.38359 


54.3 ±5.1 


2391.37270 


30.1 ± 4.3 


4169.43214 


43.8 ±1.2 


4952.39/61 


53.6 ± 3.9 


2391.38337 


32.7 ± 3.6 


4169.44290 


47.3 ± 1.4 


4952.39762 


49.6 ± 2.2 


2392.36580 


28.3 ± 2.1 


4169.45704 


56.5 ± 2.0 


4952.41400 


58.4 ± 1.8 


2392.37993 


24.5 ± 2.0 


4169.46774 


52.6 ±6.1 


4952.41401 


51.8 ±2.9 


2739.51074 


45.0 ±3.6 


4186.54793 


38.0 ±3.6 


4952.42802 


51.1 ±2.3 


2739.52029 


46.7 ± 3.8 


4186.55810 


40.1 ± 11.8 


4952.42802 


58.4 ± 1.5 


2740.44514 


49.8 ± 3.2 


4188.41925 


39.9 ± 3.0 


5202.92761 


36.5 ± 2.0 


2740.45465 


51.9 ± 3.0 


4188.43328 


42.1 ± 2.9 


5278.71693 


32.8 ± 2.0 


2741.50896 


49.7 ± 3.1 


4189.53608 


49.9 ± 1.5 


5295.96094 


35.0 ±1.6 


2741.51846 


49.4 ±3.1 


4189.53609 


42.6 ± 1.0 






* The SAAO 


science frames belonging 


; to these 4 radial velocities were made with 


very marginal observing 


conditions. The 



resulting uiu'eliable RVs are therefore not considered in the remaining analysis. 

* This is a discrepant WHT data point, taken during service observations in 2006. It is not in line with the other data taken 
at on the same night and therefore given no weight in the remaining analysis. 

* For the blue WHT science frames belonging to these two RVs, only a single arc frame was available, taken before the 
first of the two subsequent observations. 

* These two WHT observations in the Hp - Hg region had no flat fields or bias frames available. 



ive to the preliminary fit. We increased the errors by otnt =2.15 
km s^^ and ctwht =3.4 km s~^ for respectively the INT and 
WHT data. The formal NOT and HET errors were left unchanged 
since these are fibre-fed instruments which do not suffer from nor- 
mal slit-guiding errors. Also the 9 SAAO RVs were left as they 
were. Second we used the average residual from each data set to 
our best orbital fit at that point to apply offsets to the INT, WHT 
and SAAO data set|3 (respectively 0.85, -2.00 and 4.63 km s"\ 



Table lAll presents the RVs at this stage. 



predicted minus observed). Finally, we scale the errors of the en- 
tire data set multiplicatively by a factor 1.244 to obtain a value 
equal to the degrees of freedom (dof). 

Figure |2] shows the resulting radial velocity periodogram for 
PG 1018-047 in the region where we find the lowest values. 
The best solution is found around 760 days, followed by a group 
of 1 day aliases and the yearly aliases of the long period around 
250 days. An overview of the best aliases with their corresponding 
values is given in Table |4] Given the large inhomogeneity of 
the data sets, we realize that the treatment of the errors described 
above might not be perfect in all details. Therefore we have con- 
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150 




Table 5. The best orbital solution for PG 1018-047 assuming a circular 
and eccentric orbit. 



10 100 

Orbital period (d) 



Figure 2. The radial velocity periodogram for PG 1018-047, showing the 
most probable orbital periods present in the data. 



Table 4. The best orbital periods found from the radial velocity period- 
ogram for PG 1018-047. The value, mass function of the companion 
and the minimum mass of the companion obtained by assuming an sdB 
mass of 0.5 M0 are also given. 



Alias 


Period (d) 




/m (Mq) 


Ma^in (Mq) 


1 


759.80 


117 


0.17245 


0.58927 


2 


0.9987 


188 


0.00018 


0.03704 


3 


241.35 


238 


0.02602 


0.24313 


4 


267.93 


258 


0.05582 


0.34031 


5 


1.0184 


259 


0.00013 


0.03332 



sidered alternative methods to weight the data and errors as well, 
but the essential result, the long period, was unchanged throughout. 

In Table|4]we also quote the mass functions of the companion, 
calculated using 



fm — 



PKt 



2ttG 



(2) 



The minimu m mass of th e companion, assuming a typical sdB mass 
of 0.5Mq ( lHeberetal.| [l984). is also given in each case. These 
numbers will be used later to constrain the nature of the companion 
star and vice versa to show that our orbital solution is plausible 
given the properties of the secondary. 

Table[5]lists the orbital parameters for our best orbital solution 
found for PG 1018-047 and in Figure |3] the radial velocity curve 
folded on the period is shown. 

3.2.3 The secondary orbit 

Narrow metal lines of a cool companion were present in addition to 
the Balmer absorption lines from the sdB star, allowing us to obtain 
radial velocity variations for the secondary star. The blue spectra 
did not contain enough signal from the companion to be useful, 
so we focussed on the HET, NOT and red WHT/INT spectra. We 
determined the radial velocities of the secondary by means of cross 
correlation with a template spectrum ( lTonrv&Davi3ll979h . 



Circular 



Eccentric 



Porb (d) 

7 (km/s) 

e 

HJDo (d) 

KsdB (km/s) 
Kms (km/s) 



759.8±5.8 
38.2±0.5 


2453335.0±10.5 
13.0±0.8 
8.1±1.0 



755.9±5.1 
38.0±0.9 
0.246±0.052 
0±24 
2453343.0±14.7 
12.6±0.8 



Table 6. The radial velocities of the cool companion. 



Bin 



# spectra 
in bin 



Average 
phase 



RV 
(km/s) 



0.100-0.125 
0.125-0.150 
0.200 - 0.225 
0.250 - 0.275 

0.4556 

0.4585 

0.4937 

0.5370 

0.5583 
0.575 - 0.600 

0.5810 

0.6000 

0.6158 
0.625 - 0.650 
0.775 - 0.800 



2 
7 
11 

9 

HET 
HET 
HET 
HET 
HET 
2 

NOT 
HET 
HET 

2 

4 



0.1247 
0.1277 
0.2176 
0.2867 



0.5759 



0.6496 
0.7826 



35.7 ± 1.9 
32.9 ± 1.3 

29.2 ± 3.3 
29.0 ± 4.3 

36.7 ± 2.8 

34.8 ± 0.9 
39.0 ± 1.5 

41.3 ± 1.1 

41.9 ±0.7 
43.5 ± 2.2 
40.2 ± 0.6 

44.5 ± 1.1 

43.6 ± 0.6 

39.4 ± 4.9 
44.8 ± 4.9 



For the HET data we cross-correlated with 61 Cyg B (K7V), 
which was obtained as part of a different program on April 
15, 2010, using HRS in its R = 30,000 mode with a 316 
groove/millimeter cross disperser centred at 6948 A. As a result 
only seven orders from the red CCD (six, on one night) were avail- 
able for cross correlation, covering 6000 - 6760 A (orders contain- 
ing Ha or strong telluric lines were excluded). 

For the other data we constructed a template from our data 
using the average red spectrum of PG 1018-047 (the upper panel 
in Figure[TJ. Every (by eye) recognisable absorption feature of the 
secondary was fitted with a Gaussian profile, similar to the method 
described in section [3.2.11 All features possibly originating in the 
primary or of telluric origin were masked out, leaving a total of 19 
lines from the cool companion in the template (see also Table|2]l. 

Unfortunately the quality of the secondary radial velocities 
obtained from the red INT and WHT spectra were not sufficient 
to derive any trustworthy orbital period from a periodogram or to 
estimate the orbital semi-amplitude Kms for the secondary. We 
therefore used our best orbital period obtained from the primary 
RVs (Table |5]( to phase-bin the 37 medium-resolution red spectra 
before the cross-correlation routine. In total 7 out of 40 bins are 
filled. The HET and NOT spectra were left unbinned. 

In Figure |4] we have plotted phase versus radial velocity. Us- 
ing least squares fitting we estimated the projected semi-amplitude 
of the secondary to be Kms = 8.1 ± 1.0 km/s, which con- 
strains the mass ratio of the secondary to the primary star to be 
q = MMs/MsdB = 1.6 ± 0.2. The value of the fit was 13.29, 
given 15 data points minus 2 fitting parameters. 
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Figure 3. Radial velocity curve for the sdB component of PG 1018-047. We have averaged the RVs per observation run. The small lower panel 
shows the residuals. The black squares are the 7 RVs from the HET spectra, the red dot is the NOT/FIES data point, yellow squares are the 
SAAO observations, the blue left triangles and red right triangles are the respectively blue and red INT data points, and the cyan up and magenta 
down triangles coiTespond to the blue and red WHT radial velocities. 



3.2.4 Eccentricity of the orbit 

In sections 13.2.21 and 13.2.31 the analysis was based upon the as- 
sumption that PG1018-047 has a circular orbit. However, we also 
tried fitting eccentric orbits using the Markov chain Monte Carlo 
(MCMC; Gilks, Richardson & Spiegelhalter 1995) method obtain- 
ing an eccentricity of 0.246 ± 0.052; the favoured period for these 
fits decreased from 759.8 ± 5.7 days to 755.9 ±5.1 days and had 
= 152, starting from the radial velocities in Table lAll To obtain 
X^/dof = 1 we needed to scale the errors by a factor 1.153. Table 
[5]presents the orbital parameters. The argument of periapsis (cj) is 
defined as the angle in the orbital plane between the ascending node 
and the line of apsides. degrees indicates that the major axis of 
the ellipse is in the plane of sky. Note as well that HJDo is given 
as the ascending node passage minus Porb /4 to be consistent with 
our definition for the circular orbit (see eq. (O). 

Computing the F-statistic we tested whether eccentricity was 
needed in our model. Assuming a locally linear model, the variable 
^ = xiirc ~ Xecc should itsclf have a distribution with 2 dof 
(the eccentricity fit parameters), independent of the of the 114 
degrees of freedom of the eccentric orbit. 

T-1 (Xcirc Xecc)/2 t ^ D^j 
^ = Xe%c/114 = 

should then be distributed as F(2;114) under the null hypothesis 
that the orbit is circular. The chances of such a large value are very 
small, meaning we reject the null hypothesis at the 99.9% signific- 
ance level in favour of an eccentric orbit. Although formally signi- 
ficant the heterogeneous nature of our data and the limited phase 



coverage lead us to be wary of claiming a definitive detection of 
eccentricity. However, it is equally true that the orbit could be sig- 
nificantly non-circular. 

3.3 The nature of the companion 

We determined the spectral type of the companion star by minim- 
izing the residuals after subtracting different template star spectra 
from the mean red PG 1018-047 spectrum in the wavelength region 
between 6390 A and 6700 A. We masked out Ha, because the line 
contains a la rge contribution from the sdB. T his optimal subtrac- 
tion routine ( iMarsh. Robinson & Woo3ll994h is sensitive to both 
the rotational broadening Vrot sin i and the fractional contribution 
/ of the companion star to the total flux. We used template stars 
from two different origins in the spectral range FO to M9: Kurucz 
models ( Munari et al. 2005, solar composition) and 43 real star 
templates jMontes et aljl997h . All templates were prepared so that 
they had the same wavelength coverage and resolution as our nor- 
malised PG 1018-047 spectrum in the Ha region. 

Figure[5]clearly shows the difference in between the dwarf 
and (sub)giant templates. The main sequence models, which have 
the lowest x^-values, show a distinct minimum at spectral type K4 
- K6, whereas the (sub)giant models seem to converge to a spectral 
type G5. This is consistent with the set of template spectra from 
iMontes et a as well, where we find the best x^ for the K5V 

star HD 201091 (61 Cyg A), followed by a K7V and K3V star. For 
M stars the worst values are found. The best non-dwarf templates 
were both G5 stars. 
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0.2 



0.4 0.6 

Orbital Phase 



1.0 



Figure 4. The orbital fit for the secondary star The red triangles are the radial velocities from the cool companion, obtained from the interme- 
diate resolution spectroscopy. The NOT measurement is plotted with a blue dot and the HET data with black squares. The radial velocity curve 
(dashed green) is also shown. 




<^ <,% <v <f> ^ ^ ^ ^ ^ ^ ^'y ^ ^ ^ ^ 



Figure 5. Results of the optimal subtraction routine for the Kuracz tem- 
plates. The circle curve are the results for dwarf stars, whereas the tri- 
angles and squares plot the reduced versus spectral type for respect- 
ively subgiant and giant stars. 



Since PG 1018-047 was observed by 2MASS, we adopt a sim- 
ilar approach as Stark an d Wade (2003) as a complementary way to 
identify the spectral type of the companion star. We choose to focus 
on J — Ks versus B — V and J — Ks versus V — Ks and con- 
verted the Stro mgren magnitu des to the Johnson system following 
the approach of lTumen ( Il990h . The calculated colours are given in 
the lower part of table |7] 

The next step is to compare the PG 1018-047 colours to a the- 
oretical grid of colours for sdB-MS binaries while varying the frac- 
tion of light at V that arises from the companion. Defining this 



Table 7. Magnitudes and colours for PG 1018-047. The B -V colour 
is calculated using the transformation formula from Turner (1990). 



2MASS (infrared) 
(Skrutskie et al. 2006 ) 
J = 13.298 (.026) 
H = 12.980 (.027) 
Ks = 12.928 (.033) 



Stromgren (visual) 
(Wesemael et al 1992 ) 
y = 13.320 (.005) 

b-y = -0.086 (.004) 

u-b = -0.073 (.006) 

mi = 0.076 (.013) 



Calculated colours 



B-V 
J-Ks 
V-Ks 



= -0.20 (.007) 
= 0.370 (.042) 
= 0.392 (.038) 



fraction as 



f 



Vms 



F'v Fv^dB + Fvms 



we find the combined colour (e.g. for B-V) from the expression 

{B - V)sdB-t-MS = 

-2.5 log [(1 - /) ■ 10-(S-^).dB/2.5 + f . ^Q-(B-V)Ms/2.5j 

Colour indices for late-type stars are taken from ljohnsonl 1 196^. 



The typical colours for a single sdB we take from lStark and Wadg 
(2003). 

Figure |6] shows that as the fractional contribution / at V from 
the secondary increases, the B ~ V index shifts to redder values. 
On the other hand, the cooler the companion becomes, the more 
the Tfs-band dominates the J — Ks index. Swapping the B — V 
index for the V — Ks colour shows a similar trend. Zooming in on 
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Figure 6. Grid of composite colours in B — V versus J — Kg space by combining the light from a typical hot subdwarf with that of a 
population I main-sequence star assuming various fractional contributions to the total light in the V band by the late-type star. The blue dot 
marks the values for PG 1018-047. A close-up is shown in the small inset figure. The diagonal cyan line indicates the location of the population 
I main sequence. 



PG 1018-047, we find a spectral classification K3 — K6, which 
is consistent with the results from the optimal subtraction routine. 
We estimate the contribution of the secondary to be 6.1 ± 1.0 % in 
the V-band. 



4 DISCUSSION 

The role of the long period binary system PG 1018-047 in refining 
our understanding of the origin and evolution of hot subdwarf stars 
depends on whether its present orbit is circular or eccentric. 

4.1 Evolution assuming a circular orbit 

If PG 1018-047 indeed has a circular orbit, one can assume that 
tidal interaction has occurred between the sdB and the dK5 star. 
This means that from a theoretical point of view PG 1018-047 
might be an important system, as it becomes a good candidate for 
formation through the first s table Roche lobe ove rflow channel de- 
scribed by iHan et all ( |2003|) . From Figure 15 in lHanetal.1 ( l2003h 
we deduce that an sdB binary with a mid-K companion is feasible. 
However these systems are all subgiants and giants evolved from B 
stars (i.e. more massive stars at a later evolutionary state). Also, to 
have stable RLOF on the first giant branch onto a companion that is 
only 0.6-0.7 Mq and still to end up with such a long period, mass 
transfer would have to be close to conservative to avoid excessive 
angular momentum loss. If we want the mass donor to evolve in a 



Hubble time, it must have started with an initial mass of 0.8 Mq 
or greater, meaning that the initial mass of the present K star must 
have been no more than ~0.3-0.4 Mq . But then the initial mass ra- 
tio go 2 is not compatible with conservative mass transfer on the 
red giant branch, leading to a contradiction. 

Second we can consider the alter nati ve common-envel ope pre- 
scription from lNelemans et^ 1 I2OOOI) . In lNelemansI ( I2OI0I) , a pop- 
ulation of sdB stars is simulated in which the first phase of mass 
transfer c an be described by this a lternative common-envelope pre- 
scriptron jNelemans elaTlhoOlah . Interestingly, in that model a 
substantial fraction of the sdB stars with low and intermediate mass 
main sequence companions have rather large orbital periods (100- 
1000 days, see their Figure 2) and the parameters of PG 1018-047 
actually fall right in a densely populated area of the model. Thus a 
sample of these long-period sdB binaries will give another way of 
testing the outcome of the common-envelope phase. 

4.2 Evolution assuming an eccentric orbit 

IClausen&Wad3 1 (201 l|) propose that binaries similar to PG 1018- 
047 can be the remnants of original hierarchical triple systems, in 
which the inner binary has merged and evolved into the sdB star, 
and the outer (current sdB-l-MS) binary was never tidally i nteracting 
and is thus irrelevant to the production of the sdB. In the lHan et al.l 
( 2002, i2003i) merger scenario, two helium white dwarfs merge to 
make an object capable of core helium burning, and the resulting 
range of masses for the s dB star is 0.3-0.8 Mq. In the new forma- 
tion channel proposed bv lClausen & Wad3 ( 1201 ll) . a helium white 
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dwarf merges with a low-mass hydrogen-burning star whose result- 
ant total mass is ~0.6 M0 . Depending on the mixing history, this 
object has either a pre-formed helium core or is helium-enriched 
throughout, and the star experiences a greatly accelerated evolution 
to the tip of the red giant branch. Only minimal loss is required to 
remove the residual hydrogen envelope from this object at the time 
of normal degenerate helium ignition, so the expected mass range 
for the sdB star in this channel is narrow and at the "canonical" 
value ^^0.5 Mq . There is ample room inside the ~760 day orbit of 
the present-day PG 1018-047 binary system to accomodate an in- 
ner binary which underwent such evolution, and there is no cause 
for the outer orbit to have been circularised. The iClausen & Wadd 
( 1201 l|) scenario could account for stars like PG 1018-047 with long 
(and possibly eccentric) binary orbits. Stable RLOF on the other 
hand predicts perfectly circular orbits. 



Observatorio del Roque de los Muchachos of the Instituto de 
Astroffsica de Canarias. We thank PATT for their support of this 
program. We also thank the ING service scheme for obtaining 
some of the data. This paper uses observations made at the South 
African Astronomical Observatory (SAAO). A portion of the data 
reported herein was obtained with the Hobby-Eberly Telescope 
(HET), which is a joint project of the University of Texas at Austin, 
the Pennsylvania State University, Stanford University University, 
Ludwig-Maximillians-UniversitatMiinchen, and Georg-August- 
Universitat Gottingen. The HET is named in honor of its principal 
benefactors, William P. Hobby and Robert E. Eberly. This paper 
also uses observations made with the Nordic Optical Telescope, 
operated on the island of La Palma jointly by Denmark, Finland, 
Iceland, Norway, and Sweden, in the Spanish Observatorio del 
Roque de los Muchachos of the Instituto de Astroffsica de Canarias. 



5 CONCLUSIONS 

With an orbital period of 760 ± 4.7 days, PG 1018-047 is the first 
long period sdB-l-MS system for which a period has been determ- 
ined. The spectral type of the companion was found to be a K5 
dwarf star, consistent with the mass ratio Mms /MsdB = 1.6 ± 0.2 
derived from the radial velocity amplitudes of both stars. However, 
one has to note that the stated numbers are only indicative of the 
true orbital parameters, since they are sensitive to the exact uncer- 
tainties assigned to the RV data and the assumption of a circular 
versus eccentric orbit. 

At first sight PG 1018-047 is a good candidate for form- 
ation through stable Roche lobe overflow if the orbit can be 
demonstrated to be circular. The predicted number of sdB binaries 
formed through this channel amounts to the largest contribution of 
sdB binar ies according to binary population synthesis calculations 
jHan et al . 2002, 2003). At the same time, the number of known 
binaries that have been confirmed as having formed through this 
channel is small/nil compared to those formed through the com- 
mon envelope channels. Alternatively, if the common envelope is 
governed by the gamma-formalism for the common envelope, as 
used in|Nelemans (2010), there is a population of long period post- 
common envelope binaries with low-mass secondaries. Thus ob- 
serving more of these binaries can constrain the first phase of mass 
transfer. 

If the orbit turns out to be eccentric, the present binary may 
instead be the remna nt of a hierarchi c al trip le-star progenitor sys- 
tem, as outlined by Clausen & Wadd jloT lb. Further observations 
are needed to better establish the orbit of PG 1018-047. 
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APPENDIX A: RV MEASUREMENTS WITH ERROR AND 
OFFSET CORRECTION 

In table I All we report the radial velocity values for PG 1018-047 
with error and offset correction, as used for the final orbital fit res- 
ults. 

This paper has been typeset from a Tj3f/ KTgX file prepared by the 
author. 
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7. Deca et al. 

Table Al. Radial velocity measurements for PG 1018-047 with error and offset correction, as used for the final orbital fit results. 



HJD 


RV 


HJD 


RV 


HJD 


RV 


-2450000 


(km/s) 


-2450000 


(km/s) 


-2450000 


(km/s) 




INT-B 


4197.52858 


44.4± 2.9 




WHT-B 


(7= 2.15 km s~ 


1, off=-l-0.85kms-i 


4547.37872 


33.9± 2.7 


cr= 3.40 km s 


"l,off=-2.01 kms~i 


2390.51113 


35.5 ± 17.2 


4547.40665 


32.0± 2.7 


3772.55466 


33.5± 3.5 


2390.52179 


33.6 ±9.2 


4547.62144 


30.4± 2.9 


3772.68565 


32.5± 5.4 


2391.37270 


31.0 ± 4.8 


4547.64844 


33.2± 3.0 


4169.43214 


41.8±3.7 


2391.38337 


33.6 ±4.2 


4902.42320 


35.7± 4.9 


4169.44290 


45.2± 3.7 


2392.36580 


29.2 ± 3.0 


4902.44418 


37.5± 5.4 


4169.45704 


54.5± 4.0 


2392.37993 


25.3 ± 3.0 




INT-R 


4169.46774 


50.6± 7.0 


2739.51074 


45.9 ± 4.2 


a= 2.15 km s 


~1, off= -1-0.85 km s-i 


4189.53609 


40.5± 3.6 


2739.52029 


47.5 ± 3.9 


1646.47135 


32.1 ±4.9 


4189.547112 


40.9± 3.6 


2740.44514 


50.6 ± 39 


1646.47650 


21.3±4.8 


4190.40982 


44.0± 3.5 


2740.45465 


52.7 ± 3.7 


1654.46074 


25.6± 4.3 


4190.42422 


48.7± 3.5 


2741.50896 


50.5 ± 3.7 


1654.46780 


32.9± 4.3 


4191.59677 


42.1±3.6 


2741.51846 


50.3 ± 3.7 


1977.48220 


50.4± 5.7 


4588.35123 


27.7± 4.3 


2741.52795 


50.6 ± 3.7 


1978.62533 


56.0± 4.4 


4588.35648 


27.5± 4.2 


2741.53909 


54.2 ± 3.7 


1979.54044 


50. 8± 3.9 


4952.38360 


53.1± 4.4 


2741.54857 


49.7 ± 4.0 


1979.59958 


66.3± 4.0 


4952.39763 


47.6± 4.1 


2741.55807 


49.0 ± 4.8 


1979.60655 


54.8± 4.0 


4952.41340 


56.4± 3.9 


2743.42040 


51.3 ±4.9 


1982.51489 


53.0± 4.4 


4952.42802 


56.4± 3.8 


2743.43451 


61.1 ± 5.0 


1982.52185 


55. 5± 4.6 




WHT-R 


2743.52128 


53.0 ±6.6 


1982.55421 


58.7± 4.5 


a= 3.40 km s^ 


"1, off = -2.01 km s-l 


2744.35895 


55.2 ± 4.1 


1982.56118 


53.4± 4.4 


3772.55467 


24.7± 3.6 


2744.37308 


50.5 ± 3.6 


1982.60059 


58.4± 4.4 


4189.53609 


47.9± 3.8 


2744.39186 


50.2 ± 3.5 


1982.60755 


54.0± 4.5 


4189.54710 


50.1±3.7 


2744.40598 


52.3 ± 3.5 


2031.36769 


62.9± 8.6 


4190.40980 


49.9± 3.6 


2744.42529 


54.7 ± 3.3 


2031.37466 


42. 2± 8.7 


4190.42420 


51.5± 3.6 


2744.43941 


53.1 ± 3.3 


2031.38641 


50.3± 4.9 


4191.59677 


50.1±3.7 


2745.37982 


49.6 ±4.1 


2032.48243 


55.7± 3.7 


4588.35125 


32.2± 4.4 


2745.39394 


46.1 ± 3.7 


2032.49633 


54.5± 4.1 


4588.35646 


20.7± 4.3 


2745.40987 


52.2 ± 3.5 


2033.39672 


54.7± 4.6 


4952.38359 


52.3±6.2 


2745.42399 


52.7 ± 3.4 


2033.40369 


54.6± 4.8 


4952.39761 


51.6±5.2 


2746.44131 


55.4 ± 3.3 


2037.45584 


54.9± 4.4 


4952.41400 


49.8± 4.5 


2746.45544 


51.8 ± 3.3 


2037.46282 


51.8± 4.1 


4952.42802 


49.1± 4.1 


2746.47412 


56.2 ± 3.1 




SAAO 




NOT 


2746.48824 


45.7 ± 3.2 


a= km s~ 


l,off = -H4.63kms-i 


(T= km s 


-1, off = 0kms-l 




r^g g _j_ 4 2 




27 7^ 6 Q 






4186.55810 


40.9 ± 12.0 


2364.40338 


22.5± 6.6 




HET 


4188.41925 


40.8 ± 3.6 


3095.38227 


31.5± 6.6 


(T= km s 


-1, off = 0kms-l 


4188.43328 


42.9 ± 3.6 


3095.40336 


23.6± 6.6 


4441.00811 


39.68± 2.0 


4193.37192 


52.2 ± 8.1 


3098.33704 


17.5±9.8 


4469.93158 


36.00± 2.0 


4195.35718 


42.4 ± 2.9 


3101.31890 


23.4± 9.2 


4502.85094 


33.50± 2.0 


4195.37121 


37.2 ± 2.7 


3101.32956 


24.3± 8.7 


4550.71405 


26.68± 2.0 


4195.61447 


36.3 ± 3.3 


3545.21953 


56.0± 7.2 


4562.68846 


29.60± 2.0 


4195.62457 


39.0 ± 4.9 


3546.20756 


57.2± 5.8 


5202.92761 


36.47± 2.0 


4197.51456 


44.6 ± 2.8 






5278.71693 


32.75± 2.0 
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